Constraints on the Extremely High-Energy Cosmic Ray Accelerators from Classical 
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We find the general requirements, set by classical electrodynamics, to the sources of extremely 
high-energy cosmic rays (EHECRs). It is shown that the parameters of EHECR accelerators are 
strongly limited not only by the particle confinement in large-scale magnetic field or by the differ- 
ence in electric potentials (generalized Hillas criterion), but also by the synchrotron radiation, the 
electro-bremsstrahlung, or the curvature radiation of accelerated particles. Optimization of these 
requirements in terms of accelerator's size and magnetic field strength results in the ultimate lower 
limit to the overall source energy, which scales as the fifth power of attainable particle energy. Hard 
7-rays accompanying generation of EHECRs can be used as a probe for potential acceleration sites. 
We apply the results to several populations of astrophysical objects - potential EHECR sources - 
and discuss their ability to accelerate protons to 10 20 eV and beyond. A possibility to gain from 
ultrarelativistic bulk flows is emphasized, with Active Galactic Nuclei and Gamma- Ray Bursts being 
the examples. 
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I. INTRODUCTION 

The existence of EHECRs with energies above 10 20 eV 
is now a well-established fact (see, e.g. [Q and for a re- 
view). Nearly perfect isotropy of the highest energy cos- 
mic rays points to their likely extragalactic origin. On the 
other hand, the sources must be located within 100 Mpc: 
the propagation distance for EHECRs of energy greater 
than 5 x 10 19 eV is limited due to their interaction with 
the microwave background radiation || ||], unless there 
are some yet unknown type(s) of particles that do not 
undergo such an interaction. 

The presence of non-thermal particles is a usual fea- 
ture of any plasma system. However, finding a way to 
accelerate these particles to 10 20 eV is a serious theo- 
retical challenge. Even if an electrodynamical system 
used to accelerate charged particles is sufficiently long- 
lived, the attainable energy is limited since the parti- 
cles can only stay within the acceleration region until 
their gyroradius exceeds the size of the region. (This is 
known as Hillas criterion ||.) Various types of model- 
and source-dependent radiative and collisional losses de- 
crease the limiting particle energy furthermore. But even 
if one finds a way to circumvent the latter obstacle, there 
are radiative losses caused by the same electromagnetic 
field that is used to accelerate particles. For example, 
Greisen || considered the synchrotron losses and showed 
that for very small-sized accelerators, i.e., for those hav- 
ing very strong magnetic field, the synchrotron radiation 
drains energy so fast that even a particle accelerated at 
the largest possible rate for the magnetic-field dominated 
environments, qBc {q is the particle charge, B the mag- 
netic field strength, c the velocity of light), is not able 
to achieve the energy given by the Hillas criterion. How- 
ever, Greisen's limit might seem to be not a general one, 
as in some objects accelerated particles move along the 



magnetic field lines (neutron stars are a perfect example 
of that), without suffering from the synchrotron losses. 

In this paper, we relax the requirements to the accel- 
eration sites as much as possible, assuming that inho- 
mogeneity of electric and magnetic fields is small enough 
not to cause particle diffusion out of the acceleration site 
and taking into account only the radiation mechanism 
which gives the lowest energy loss rate. In fact, even if 
the synchrotron losses (or the electro-bremsstrahlung in 
the case of electric-field dominated environments) tend to 
zero, i.e., for particles moving along smooth field lines, 
there is still curvature radiation, which then limits the 
maximal particle energy. Furthermore, we generalize the 
analysis introducing acceleration efficiency parameter rj, 
so that the acceleration rate is given by 



e = r/qBc , 



(1) 



where rjB = E e g is the projection of an electric field on 
the particle trajectory, effectively averaged as the parti- 
cle moves along this trajectory. Such a representation is 
useful when, e.g., diffusive shock acceleration scenario is 
brought into consideration. 

We derive the ultimate lower limit for the total elec- 
tromagnetic energy stored in an acceleration region and 
find the size and the field strength that are optimal with 
respect to minimization of this energy. Then we consider 
a possible gain from relativistic motion of an accelerator 
as a whole and discuss what kind of additional informa- 
tion one gets from analysis of hard 7-photons generated 
as a byproduct of cosmic ray acceleration. Finally, we 
make comparison of theoretical constraints with the pa- 
rameters of known astrophysical objects. 
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II. ELECTRODYNAMICAL LIMITATIONS 

Let us consider a region of space of a size R where par- 
ticle acceleration occurs. A natural assumption is that 
there is nothing like linear accelerator geometry in this 
region, but instead the field lines of both electric and 
magnetic fields are curved with the curvature radius be- 
ing of the order of R or less. 

If the magnetic field is much stronger than the electric 
field, B 3> E, then a test particle is accelerated follow- 
ing a magnetic field line in its trajectory. If E > B but 

(E ■ B) ~ EB, then particle's trajectory depends on the 
radiative loss rate. High loss rate will force the particle 
to follow a field line (of an electric field in this case) since 
a transverse to the field line component of particle's mo- 
mentum will quickly decay due to an increase in radiative 
losses. In the opposite case, the field geometry will have 
little effect on particle's trajectory which therefore will 
make an angle of ~ 45° with field lines. The remaining 
case, where both Lorentz invariants (E-B) and \E 2 — B 2 \ 
are much smaller than E 2 + B 2 , can be reduced to one 
of the already mentioned situations by the appropriate 
Lorentz transformation to the reference frame in which 
the mean (after averaging in space and time) Poynting 
flux density becomes zero. 

In what follows we assume that accelerator's lifetime 
is at least as large as the acceleration timescale and use 
the most optimistic premise that the strongest of mag- 
netic and electric fields is regular on scales <C R and 
quasistatic. (Alternating field direction implies that par- 
ticles escape more readily and that the time-averaged 
acceleration rate is smaller, while the radiative losses are 
roughly the same.) Then, the effect of small-scale mag- 
netic field can be neglected if E > B, and in the case 
B > E the factor i] < 1 takes into account all necessary 
averaging for a small-scale and/or variable electric field, 
so that E e g = rjB. 

The particle can be accelerated up to the terminal en- 
ergy which is the smallest of either the work done by 
accelerating force in a time it takes for a particle to es- 
cape from the region, or the energy at which radiative 
losses balance the acceleration. The balance condition 
requires that the radiative loss rate e ra( j is equal to the 
energy gain rate, i.e. 



2 4 q 2 
£rad = 77 T^ c = m Bc 



(2) 



3 ' R 2 

in the case of curvature-radiation dominated losses 
(the particle is assumed to move along an arc of radius 
R), or 



2 2 ( q 2 N 1 



£rad = 7/7 I 2 

3 \ mc 



c (B 2 ± + El) = rjqBc (3) 



in the case of synchrotron- or electro-bremsstrahlung 
dominated losses. Here 7 is the Lorentz factor of the 
particle, m its mass, B± and E± are the r.m.s. field 
components perpendicular to the particle momentum. 



The comparison of Eqs. (g) and (||) shows that the 
curvature radiation gives a more favorable (smaller) esti- 
mate for the radiative losses as long as the particle energy, 
e = jmc 2 , satisfies the condition 



e < qRJBi 



E 2 ± . 



(4) 



If the above condition were violated, it would mean one 
of the following. Either the particle is no longer bound in 
the acceleration region, i.e., its gyro-radius rjj = sj {qB±) 
is larger than R (when rn ^ R, the very idea of cur- 
vature radiation becomes meaningless) and acceleration 
terminates upon the particle escape at e < qRB±, or, 
if B±_ -C E±, the particle is unbound from the very 
beginning and its energy gain is limited by the differ- 
ence of electric potentials across the acceleration region, 
e < qRE c ff = rjqRB. This provides a generalized version 
of the Hillas criterion. 

Hereafter we neglect numerical factors of ~ 1 assuming 
B = B±_ and E = E± = E e s = rjB. In general, one 
should bear in mind that for an inhomogeneous small- 
scale field the r.m.s. quantity E 2 must be replaced with 
(E 2 ) = E 2 /5, where 5 < 1 is a duty-factor. (The same 
is true for magnetic field, although the situations where 
E > B and the magnetic field has the duty-factor 5<1 
seem quite unrealistic.) We assume below 6=1, aiming 
to the most efficient acceleration. Thus, the condition 
given by Eq. (Q) cannot be violated, as it satisfies the 
generalized Hillas criterion, 



£max < qRmax{B,E}. 



(5) 



Since we are interested in the most efficient acceler- 
ation (i.e., the highest acceleration rate at smallest ra- 
diative losses), it is Eq. (||) that will be used below to 
put the lower limit on the radiative loss rate. Then, we 
conclude that the maximal energy of a particle is either 
determined by the balance condition given by Eq. (0), 



3rjBR 2 

2g 



mc 



(6) 



or is limited by the Eq. (^J), whichever gives the smaller 
value. 

The overall energy of an electromagnetic field in a 
spherical region of radius R, capable of accelerating par- 
ticles up to e max , is Wem = R 3 (B 2 + E 2 )/Q. According 
to Eqs. (||) and respectively, it is limited as follows: 



2 q 2 /£max\ 8 l + V 2 



Wcm > 27 R V mc 2 ) 



and 



W cm > — 



R (e 



(7a) 



(7b) 



Our interest is to lower the energy requirements taking 
R as large as possible in Eq. (|7|a) and as small as possible 
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FIG. 1: The energy requirements for protons of energy 
10 20 eV. Eq. (fjja), dotted lines (for acceleration efficiencies 
r] — 1 and r\ — 10~ 2 ), and Eq. (j?]b), dashed line, as functions 
of the size of an acceleration region. 



m Eq. (0b). Doing this, we derive the optimal size 



r?(°pt) ~ _ 
" 3 



2^1 



V 2 <2 2 £max 



(mc 2 ) 



(8) 



such that both conditions become quantitatively the 
same. This finally gives the optimal (the minimum pos- 
sible) estimate for the amount of electromagnetic energy 
stored in the acceleration region: 



I4/(°pt) 



1 a/1 + rf 
9 



(9) 



V (mc 2 ) 

The corresponding optimal strength of the magnetic field 
is 

BW^JL-fci!, (io) 
2 1 + ry 2 g3 £ 2 nax 

and E(° p Q = T]B^ opt \ 

The energy requirements given by Eqs. (0) as functions 
of accelerator size are presented in Fig. 1. The required 
accelerator energy increases as the region size deviates 
from the optimal one. A smaller region would accelerate 
particles to the energy £ max only if the energy budget 
is mcr eased by R^/R times; in a larger region, the 
energy budget has to be increased by R/R^° pt ^ times. 

To have an idea of typical numbers, let us make the 
estimates for a proton of energy e max = 10 20 eV. In this 
case, the total electromagnetic energy in the acceleration 
region should exceed Wem^ — 3 x 10 51 erg; the optimal 
size i?( opt ) is of the order of 10 17 cm and the optimal mag- 
netic field is of ~ 3 G. When making these estimates we 
assumed rj = 1 to lower the energy requirements. Taking 
X] still larger does not reduce them too much, by a factor 
of 2 at most. 

All the results obtained arc valid also for an accelera- 
tor that moves as a whole with the Lorentz factor r>l, 
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FIG. 2: The energy requirements for protons of energy 
10 20 eV. Eq. (Qa), dotted lines (for accelerator at rest and 
for bulk Lorentz factors F = 10 and T = 300), and Eq. (Qb), 
dashed lines, as functions of the comoving size of an accelera- 
tion region. For AGN inner jets, the upper zone corresponds 
to hadronic models of high-energy 7-ray emission, and lower 
zone - to leptonic models. 



if the quantities are measured in the comoving frame. 
However, it is more convenient to present the results in a 
form where Wem and e max are measured in the labora- 
tory frame, while R^ opt \ B^ opt ^ and £;( opt ) are measured 
in the comoving frame. Then 



^>/(opt) 



£/(opt) 



1 VT+ 



9F 4 77 


(mc 2 ) 4 


2 + 


V 2 9 2 £max 


3F 3 7] 


(mc 2 ) 


3T 2 r, 


(mc 2 ) 4 



2 1 



? '-max 



(11a) 



(lib) 



(11c) 



where the primed quantities are measured in the comov- 
ing frame. A comparison of energy requirements given by 
Eqs. (Q), generalized for the case V » 1, for accelerators 
with different bulk Lorentz factors is presented in Fig. 2. 

Apparently, the acceleration of EHECRs in ultrarela- 
tivistic bulk flows has an advantage of bringing the opti- 
mal source size to a more "comfortable" (for short-lived 
phenomena like Gamma-Ray Bursts) range and reduc- 
ing the energy requirements to the source. However, for 
wind-like relativistic flows the actual requirements are 
geometry-dependent. Indeed, for the causality reasons, 
the acceleration region does not occupy the whole sphere 
of radius R, but rather extends to a distance R' = R/T 
transverse to the radius and to a distance R' /T along it, 
so that the total energy within the volume of radius R 
is of ~ r 4 Wo m for a wind with a broad beam pattern. 
Since the energy stored in the acceleration region, W em , 
is multiplied by T 4 in this case, one can only gain from 
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a more favorable ratio R'/R'^ (cf. Eq. (0a)). A jet 
geometry is more favorable: with the jet opening angle 
~ 1/r the energy requirements can be reduced by a fac- 
tor as large as T 2 , with respect to the non-relativistic 
accelerator. To implement the advantage of relativistic 
bulk motion to the full extent, i.e., to get a gain by a 
factor of r 4 , one needs a transient jet source which lasts 
for a time ~ R/{T 2 c). 



III. CONSTRAINTS DUE TO RADIATION 
FROM ACCELERATION SITES 



Before moving to an overview of different source pop- 
ulations, let us discuss some implications regarding the 
acceleration timescale and the electromagnetic radiation 
accompanying the acceleration process. 

When r] > 1, a particle is accelerated up to e max in 
a time i acc < R/c, which is less than or comparable to 
the dynamical timescale of an accelerator. However, at 
low efficiency ()<1, the acceleration may take a longer 
time, i.e. rj~ 1 {R/c), if accelerator's size is larger than 
i?(°Pt), or r]- 1 (R/R ( - op ^)(R/c) otherwise. Relativistic 
motion of an accelerator usually implies relativistic ex- 
pansion. Consequently, the lifetime is of order R/c, so 
that small r\ strongly disfavors such objects unless their 
size is < riR(° pt \ that would allow a particle to attain 
the energy e max . 

The acceleration of particles goes side by side with the 
emission of hard electromagnetic radiation. The lower 
limit on the energy lost to this radiation is given by the 
product of t acc and the radiative loss rate (Eq. (|^) , where 
it is convenient to substitute 7 3 from Eq. (||)). However, 
if iacc < R/c one must integrate the loss rate over a time 
of order R/c, since a particle cannot leave the acceler- 
ation region in a shorter time. As a result of particle's 
acceleration, the following fraction of its terminal energy 

£max is lost for 7-radiation 



fl(opt) j R if i? > R(°pV 

} 1 if i?(°pt) > r > rjR(opt) 

rjR^/y/l + ^R if R<r)R(° pt \ 

(12) 

Therefore, acceleration in small objects (R < i?( opt )) is 
always limited by radiative losses. 

The spectral energy distribution of the accompanying 
radiation has a maximum (in the source frame) around 



7 3 c 



R aZ 2 ' 



(13) 



provided the accelerator size is less than R {o ^; a is the 
fine-structure constant and Z the ion charge in units of 
the elementary charge. The energy of curvature radiation 
photons given by Eq. ( |l3| ) is related to the energy of 
synchrotron photons as e 7 = (rn/R)e^. 

The ratio of the luminosity in accompanying radiation 
to the luminosity in accelerated particles, L 7 /L acc , as 



well as the energy of individual photons, is reduced if the 
acceleration time is limited to R/c, like in the sources 
with relativistic bulk motion (particle diffusion may lead 
to analogous limit), but this is at expense of reducing 
£max- Since the radiative loss rate (Eq. (j^)) and the 
energy of individual photons (Eq. (J13|)) scale as 7 4 and 
7 3 , while particle's Lorentz factor increases linearly with 
time, the corresponding factors for L 7 /L 



7/ -^acc 5 °7 



and £„ 



are: (rjR^ / R) 4 , {r/R^/R) 3 , and (rjR^/R) respec- 
tively, in the case where i? (opt) > R > ?7i? (opt) . Note that 
if the maximal energy were limited by the synchrotron 
rather than curvature losses, then the acceleration up to 
£max would not be possible at all in the objects smaller 
than R(°pV if r) < 1, as noted in @. 

The check whether the sky-averaged luminosity in ac- 
companying electromagnetic emission exceeds the diffuse 
background is a powerful selection tool for possible ac- 
celeration sites. For example, the acceleration of protons 
with high efficiency, 77 ~ 1, results in 7-radiation in en- 
ergy range > 100 GeV (Eq. (|l3|)). Such photons un- 
dergo reprocessing into 10 MeV - 100 GeV energy range 
through interaction with cosmic infrared and microwave 
backgrounds. The observed spectral energy distribution 
of diffuse 7-ray background, ~ 10 -9 erg/ (cm 2 s sr) _1 [||, 
is virtually constant in the 10 MeV - 100 GeV spectral 
domain, and the observed flux of 10 20 eV cosmic rays 
is - 2 - 5 x 10~ 12 erg/(cm 2 s sr)" 1 §, i.e., 200 - 500 
times less. But the GZK effect limits the distance to 
which the sources of 10 20 eV cosmic rays can be ob- 
served to ~ 30 Mpc, while the related 7-radiation is vis- 
ible up to a redshift z ~ 1, i.e. the distance ~ 3 Gpc. 
If the accelerators are uniformly distributed in the Uni- 
verse and no other contributions to the 7-ray background 
are taken into account, this implies L 7 /L acc < 5 and 
hence R/(rjR^) > 1/5. Since the spectral energy dis- 
tribution of the 7-ray background is nearly constant in 
10 MeV - 100 GeV range, the same limit applies for iron 
nuclei, although their individual 7-luminosity is concen- 
trated in the lower-energy spectral domain > 10 GeV 
(see Eq. (|l|), which gives scaling e 1 oc A/Z 2 ). 

The accompanying radiation may also help in identi- 
fying acceleration sites. The electrons are more suitable 
for this purpose as they can achieve radiative-loss limited 
regime of acceleration in an object with much smaller en- 
ergy content. For example, there is an excess emission in 
100 MeV range from the Crab plerion, which was inter- 
preted as the electron synchrotron radiation Q . Interest- 
ingly, the size of the plerion, ~ 10 17 cm, and its estimated 
magnetic field, ~ 10 -4 G, are very close to the optimal 
values for electrons accelerated at the efficiency 77 = 1. 
The electrons then attain the energy ~ 3 x 10 15 eV and 
the peak of their synchrotron spectrum is at ~ 100 MeV. 

So far, we did not take into account the proton interac- 
tion with smallest-scale electromagnetic field, which re- 
sults in photomeson reactions (at wavelengths < 0.04 cm 
or photon energies > 2 x 10~ 3 eV for 10 20 eV pro- 
ton). Such an interaction, with the cross-section a pl ~ 
2 x 10~ 28 cm 2 and inelasticity / ~ 0.5, may in some cases 
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lead to unacceptably high energy losses for the acceler- 
ated particles. For illustration, we make the estimate for 
10 20 eV protons and for the most favorable case 77 = 1. 
In order to allow accelerated protons to escape without 
significant losses, the luminosity of the source at the max- 
imum of integral photon spectrum must be 



Rc 

L < IT— £ ■ 



fa 



1.5 x 10 



44 



in 



R 



10 17 < 



( Tiv) erg/s 



(14) 



where e is the location of the maximum. 



IV. OVERVIEW OF POSSIBLE 
ACCELERATION SITES 

Neutron stars. Because of their size, which is far 
too small compared to R^° pt \ the performance of neu- 
tron stars as cosmic ray accelerators is limited be the 
curvature losses (see, e.g., @). For 10 20 eV EHECRs, 
the energy requirements to neutron star magnetosphcrcs 
are unphysical: ~ 10 62 erg for protons and ~ 10 51 erg 
for iron nuclei. Acceleration of heavier nuclei, with mass 
A ~ 150 and charge Z ~ 50, to slightly smaller en- 
ergy 5 x 10 19 eV is marginally possible for rapidly ro- 
tating (i.e. rj ~ 1) neutron stars with magnetic field 
B ~ 10 15 G (so-called magnetars). However, the ques- 
tion whether such nuclei can be efficiently transported to 
neutron star surface and ejected from it remains open. In 
addition, the luminosity in accompanying 10 — 100 TeV 
curvature 7-radiation (which is likely to be reprocessed 
through electromagnetic cascade into 1 — 100 MeV range) 
would be 5 to 6 orders of magnitude higher than in the 
produced EHECRs. Particles accelerated near the light 
cylinder produce less accompanying radiation, but the 
energy requirements rise in this case. Therefore, magne- 
tospheres of neutron stars have to be excluded from the 
list of EHECR sources. 

A way out from the problem of catastrophic overpro- 
duction of the accompanying radiation is acceleration 
of particles within ultrarelativistic pulsar winds. How- 
ever, the energy constraints, which now have to be re- 
formulated in terms of required Poynting-flux luminosity 
L cm ~ T 4 W em c/R for a wind with a broad beam pattern, 
remain extremely tough. The best case (R' — i?'( opt ) and 
rj = 1) for 10 20 eV protons corresponds to T < 600 and 
L em ~ 10 45 r 2 erg/s. With iron nuclei the figures be- 
come r < 60 and L cm ~ 1.5 ■ 10 42 r 2 erg/s. In any case, 
the required luminosity is far beyond any of the observed 
values. 

Black holes. Rapidly rotating black holes embedded 
in a magnetic field can generate electric field and accel- 
erate particles (e.g., ]h| [n]). We make an (optimistic) 
estimate assuming that the magnetic field is roughly in 
equipartition with accreted plasma, an accretion rate cor- 
responds to the Eddington regime, and acceleration oc- 
curs within 3R g , where R g ~ 3 x 10 5 cm M/Mq is 
the Schwarzschild radius of the black hole of mass M. 



R °> 81^ 



<7t- 



The energy W em ~ nm v c R g /aT is available in the ac- 
celeration region (m p is the proton mass, and <7t — 
6.6 x 10 -24 cm 2 the Thomson cross-section). This is to be 

compared with the required energy (i?( opt )/3 Rg)Wem , 
which yields 

)l/3 
(15) 

~3 x 10 15 Z 2 / 3 A- 8 / 3 cm. 

If the accelerated particles are protons, then only super- 
massive black holes with M > 1O 1O M0 can meet the 
above requirement. On the other hand, a lighter hole of 
M > 2 x 10 6 M Q can accelerate iron nuclei, and, since the 
size of such a hole is close to the optimal one or exceeds 
it (R g > i?(°P*) ~ 8 x 10 12 cm), this would not be at odds 
with the observed 7-ray background. 

Active Galactic Nuclei (AGNs). This phe- 
nomenon is due to ejection of a relativistic jet by a mas- 
sive black hole interacting with an accretion disk around 
it. Inner jets in AGNs have typically Lorentz factors 
between 5 and 10. The estimate of the magnetic field 
depend on the radiation model. For leptonic models 
which imply high-energy 7-ray production through in- 
verse Compton scattering, the magnetic field is of ~ 0.1 G 
(e.g., iQ). However, in hadronic models, in particular in 
the extreme proton-synchrotron model, the field strength 
exceeds 100 G O]. Furthermore, the size of accelera- 
tion region, 10 10 16 cm as inferred from variability 
timescale, is close to the optimal value for 10 20 eV pro- 
tons. This eliminates problems with 7-ray background 
and minimizes energy requirements that constrain the 
Poynting-flux luminosity of the source in the case of a 
relativistic outflow. This luminosity (rW em c/ff) must 
be higher than ~ 10 45 erg/s (or a factor of (6T) 2 more if 
the opening angle of the jet, 0, is larger than T -1 ), which 
is reasonable for a typical AGN. 

Large accelerators. There are three kinds of very 
large objects that were proposed to produce EHECRs: 
knots (or hot spots) in large-scale AGN jets, radio lobes 
(shocked regions), and galaxy clusters (e.g., || [l4|, p"5|). 
Despite their sizes are far from the optimum: ~ 1 kpc, 
~ 100 kpc, and ~ 1 Mpc respectively, the energy of 
a magnetic field is enough to bind 10 20 eV protons. 
The limits on these objects are better formulated in 
terms of the requirements to the acceleration efficiency 
77, and they are likely to arise from particle diffusion. 
In all cases there is about 1000-fold excess of the en- 
ergy stored in magnetic field with respect to that shown 
in Fig. 1. Therefore, if particles escape on the Bohm- 
diffusion timescale ~ 3 R 2 / (rnc), a rather high acceler- 
ation efficiency 77 > 0.03 is necessary. Even if the diffu- 
sion is strongly suppressed (e.g., special configuration of 
magnetic field lines), there are limits on rj. It must be 
larger than ~ (0.3 — 1) x 10~ 2 in radio lobes and galaxy 
clusters, otherwise protons are accelerated so slowly that 
they start to lose energy interacting with the cosmic mi- 
crowave background, which eventually stops acceleration. 
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In hot spots, with their stronger magnetic field up to 
B ~ 10~ 3 G, r) can be as small as few xl0~ 5 , in which 
case their size becomes slightly smaller than the optimal 
one leading to a further increase in energy requirements. 

Gamma-Ray Bursts (GRBs). These explosion- 
like events have been argued to be suitable places for 
acceleration of EHECRs (e.g., jl6], [l?]]), owing to their 
compactness (variability timescale less than a second), 
high apparent (47r) luminosity of the order of 10 52 erg/s, 
and ultrarelativistic bulk flows with the Lorentz factors 
r ~ 100 — 1000, which are the key ingredients of any 
GRB model (see, e.g., Q for a review). Acceleration 
may take place either in internal shocks formed within 
the outflow (at a distance 10 10 — 10 13 cm from the cen- 
tral engine) by colliding shells of different Lorentz factors 
or in the external shock (at a distance 10 16 — 10 17 cm) 
at the interface between ejected material and interstellar 
gas. The distance from which the accelerated protons can 
escape is, however, limited at least by their interaction 
with the intense thermal photon field left after adiabatic 
expansion of relativistic fireball. The escape distance is 
(in the shock comoving frame) 



0.08 



fa pi L _ 2 x 10 17 cm 
7iT 3 cT ~ 



3 MeV 



T 3 10 52 erg/s T 



(16) 



where L is the GRB 4 7r-luminosity and T the temper- 
ature at the base of the fireball, T ~ 3 MeV. The 
Eqs. (fj"l|b) and ( |l6| ) scale with the same power of T, 
so that the ratio between R' csc and i?'( P*) depends prac- 
tically only on GRB luminosity, and the escape distance 
appears to be larger than the optimal size, unless a 
GRB has unusually low luminosity. Therefore, radiative 
losses are negligible and the maximal particle energy is 
e max = rjqB' c(R' / c) . With the magnetic field related to 
the GRB 4 7r-luminosity via the parameter e m , which is 
the fraction of energy carried in the form of magnetic 
field, this gives 



= IV 

max 



rjc 



1/2 



2L 2.5xl0 23 eV 



r 



(17) 



Note that the maximal energy e max decreases with in- 
creasing bulk Lorentz factor T. This can be understood 
following the same reasoning as at the end of the Sect. II: 
changing the Lorentz factor allows one to gain only from 
a more favorable ratio i?'( opt )/i?', but for GRBs the es- 
cape distance is larger than the optimal size, so that the 
increase in L makes things worse. 

Since explanation of some GRB properties, such as 
rapid variability and the presence of emission in the 
10 GeV domain 113], requires L > 100, we conclude 
that GRBs might be capable of accelerating protons well 
above 10 20 eV. However, to match the observed cosmic 
ray flux at least 50 per cent of the observed GRB energy 
has to be converted into EHECRs. 



V. CONCLUSION 

Radiation losses, either due to synchrotron or curva- 
ture radiation, play increasing role in the energy bal- 
ance of accelerated particles with decreasing size (and 
accordingly increasing magnetic field) of the accelerator. 
In sufficiently compact objects the radiative losses limit 
maximal particle energy to the values that are smaller 
than those determined from the Hillas criterion. Thus, 
there are optimal size and magnetic field strength that 
minimize the requirements to the internal energy of the 
accelerator. For 10 20 eV protons and acceleration effi- 
ciency 77 = 1, the optimal size is i?( opt ) ~ 10 17 cm and 
the optimal magnetic field is B^ ~ 3 G. This corre- 
sponds to the total electromagnetic energy in the accel- 
eration region W&n — 3 x 10 51 erg. Relativistic bulk 
motion reduces the energy requirements and the opti- 
mal comoving size, while increasing the optimal comov- 
ing magnetic field. The Lorentz factor T — 10 leads to 
#'(°pt) ~ iqU cm an( j B '(opt) ^ 3 00 q corresp onding to 

the total electromagnetic energy W e n? — 3 x 10 47 erg 
in the region of size i?'( opt ) . 

Efficient acceleration in optimally-sized objects leads 
to significant radiative losses, which in this case occur 
in a regime transitional between the synchrotron and 
the curvature radiation. The emission is concentrated 
at photon energies of ~ 100 — 300 GeV (protons) and 
~ 50—150 MeV (electrons) times the bulk Lorentz factor. 
This suggests the use of new generation of air-Cherenkov 
telescopes and GLAST, respectively, for a search for ef- 
ficient accelerators. 

The derived constraints seem to firmly exclude any 
known population of neutron stars as acceleration sites 
of EHECRs. Black holes also appear to be not able to 
accelerate protons to 10 20 eV except for the supermassive 
ones with masses M > 1O 1O M , although acceleration of 
heavy nuclei by lighter holes is still a possibility. Other 
astrophysical objects considered above are able, in prin- 
ciple, to accelerate EHECRs, although the required effi- 
ciency r) is always comparable to unity. In the standard 
treatment of the problem, such an efficiency is difficult 



to achieve even in ultrarelativistic shocks 20 . In ad- 



dition, energy losses due to photomeson reactions may 
be a serious problem for radiation-rich environments of 
supermassive black holes and for GRB sources. 

The results obtained in this paper are the unavoidable 
limits, set by the laws of classical electrodynamics. Tak- 
ing into account the effects of diffusion, the energy losses 
via photomeson process, the energy content of particles 
and waves in cosmic plasma, etc., should lead to still 
stronger limits. However, they cannot, of course, con- 
strain the top-down models for the origin of EHECRs 
(e.g., pl[ ), which make use of quantum effects and non- 
electromagnetic interactions. 

In view of the development of advanced experiments 
which aim at detection of cosmic rays (CRs) above 
10 21 eV, it is worth noting what is the highest energy of 
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EHECRs still consistent with the acceleration scenario of 
their origin. Plotting diagrams similar to Figs. 1 and 2, 
but for different values of e max , one can see that in GRBs 
and AGN inner jets protons can in principle be acceler- 
ated up to 10 21 eV. At such energies, acceleration is likely 
to be radiative-loss limited in both cases. Above 10 21 eV 
GRBs fail because of their insufficient duration and AGN 
inner jets - because of insufficient Poynting-flux luminos- 
ity. Hot spots, radio lobes and galaxy clusters can still 
work to (3 — 5) x 10 21 eV under very speculative assump- 
tion that the magnetic field is ordered on all scales and 
the acceleration efficiency is rj ~ 1. In this 
celcration is escape-limited. At the energies > fO 22 eV 



the CR primaries have to be heavy nuclei. In all the 
sources listed above the heavy nuclei are accelerated in 
the escape-limited regime, so that the attainable energy 
is roughly Z times more than for protons. However, the 
nuclei of such energy are fragmented through interaction 
with the microwave background photons after traveling 
a distance of less than I Mpc, that means they must be 
produced within the local group of galaxies, and GRBs 
is the only possibility to do this. Although the nuclei 
are easily fragmented in radiation-reach environments of 
GRBs, we have to conclude that formally the primaries 
with energy up to (2 — 3) x 10 22 eV can be produced 
within the framework of acceleration scenario. 
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